R Notebook

library(raster)
library(sf)
library(tmap) # for making the nice plot
library(osmdata)
library(OpenStreetMap)
library(tidyverse)
library(terra)
library(spDataLarge)
library(stars)
library(ggplot2)
library(dplyr)
library(tmaptools)

ontario_bbox <- c(left = -95.1537, right = -74.3202, bottom = 41.6766, top = 56.8566)

# Query OpenStreetMap for administrative boundaries of Toronto
ontario_border <- opq(bbox = ontario_bbox, timeout=120) |>
  add_osm_feature(key = "boundary", value = "administrative") |> 
  add_osm_feature(key = "ISO3166-2", value = "CA-ON") |>
  add_osm_feature(key = "admin_level", value = "4") |> 
  add_osm_feature(key = "type", value = "boundary") |> 
  add_osm_feature(key = "wikidata", value = "Q1904") |> 
  add_osm_feature(key = "name:en", value = "Ontario") |> 
  osmdata_sf() |>
  (\(x) x$osm_multipolygons)() |>
  select(osm_id,admin_level,border_type,name,boundary)

ontario_border = st_transform(ontario_border, crs = 4326)
ontario_basemap = read_osm(ontario_border, type = "esri-topo")


ontario_border = st_transform(ontario_border, st_crs(ontario_basemap))
ontario_cropped <- st_crop(ontario_basemap, ontario_border)


tm_shape(ontario_cropped) +
  tm_rgb()

brampton_bbox <- c(left = -79.8587, right = -79.6633, bottom = 43.6319, top = 43.7559)

# Query OpenStreetMap for administrative boundaries of Toronto
brampton_border <- opq(bbox = brampton_bbox, timeout=120) |>
  add_osm_feature(key = "boundary", value = "administrative") |> 
  add_osm_feature(key = "admin_level", value = "8") |> 
  add_osm_feature(key = "type", value = "boundary") |> 
  add_osm_feature(key = "wikidata", value = "Q44198") |> 
  add_osm_feature(key = "name", value = "Brampton") |> 
  osmdata_sf() |>
  (\(x) x$osm_multipolygons)() 

brampton_border = st_transform(brampton_border, crs = 4326)
brampton_basemap = read_osm(brampton_border, type = "esri-topo")


brampton_border = st_transform(brampton_border, st_crs(brampton_basemap))
brampton_cropped <- st_crop(brampton_basemap, brampton_border)

tm_shape(brampton_cropped) +
  tm_rgb()

map_base <- tm_shape(brampton_cropped) +
  tm_rgb() +
  tm_layout(title = "Brampton Community Projects Base Layer")

print(map_base)


map_with_first <- map_base +
  tm_shape(brampton_first) +
  tm_dots(col = "red", size = 0.1, popup.vars = c("id", "name"), legend.show = TRUE, title = "2014-17 Projects")


print(map_with_first)


map_with_second <- map_with_first +
  tm_shape(brampton_second) +
  tm_dots(col = "blue", size = 0.05, popup.vars = c("id", "name"), legend.show = TRUE, title = "2017-20 Projects")

print(map_with_second)

map_with_third <- map_with_second +
  tm_shape(brampton_third) +
  tm_dots(col = "green", size = 0.025, popup.vars = c("id", "name"), legend.show = TRUE, title = "2020-24 Projects") +
  tm_layout(title = "Brampton Community Projects 2014-24",
            legend.outside = TRUE, 
            legend.outside.position = "right",  # Position the legend outside the plot to the right
            legend.title.size = 1.2,  
            legend.text.size = 0.8,  
            legend.bg.color = "white",  
            legend.bg.alpha = 0.5)


print(map_with_third)

tmap_arrange(healthcare_brampton_first_plot, healthcare_brampton_second_plot, healthcare_brampton_third_plot, ncol = 3)

Error: The shape healthcare_brampton_second only contains empty units.